Experimental implementation of assisted quantum adiabatic passage in a single spin 
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Quantum adiabatic passages can be greatly accelerated by a suitable control field, called a counter- 
diabatic field, which varies during the scan through resonance. Here, we implement this technique 
on the electron spin of a single nitrogen-vacancy center in diamond. We demonstrate two versions 
of this scheme. The first follows closely the procedure originally proposed by Demirplak and Rice 
(J. Phys. Chem. A 107, 9937 (2003)). In the second scheme, we use a control field whose amplitude 
is constant, but its phase varies with time. This version, which we call the rapid-scan approach, 
allows an even faster passage through resonance and therefore makes it applicable also for systems 
with shorter decoherence times. 

PACS numbers: 03.67, 33.35, 76.70 



Introduction. - Controlling quantum systems with high 
fidelity is an essential prerequisite in various fields, such 
as coherent control of atomic and molecular systems [l| 
and quantum information processing 0, 0] . The strate- 
gies that have been developed for this purpose include 
the adiabatic passage technique, which leads the quan- 
tum system along a specific pathway in such a way that 
the system always remains in its ground state. One of 
the attractive properties of this technique is that the re- 
sulting evolution is robust with respect to some experi- 
mental imperfections Adiabatic passage also is the 
central part of the adiabatic model of quantum compu- 
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tation [5|, |6|J , which has been shown to be equivalent to 
the more common network model. In all these cases, it is 
essential that the scan duration of the adiabatic passage 
is short and the fidelity as high as possible. 

The quantum adiabatic theorem guarantees that the 
system remains approximately in its ground state if the 
evolution is sufficiently slow However, for all prac- 

tical applications, the optimal implementation is reached 
when the scan time remains short, e.g. compared to the 
decoherence time. This conflict can be relieved by an ad- 
ditional control field, the so-called the counter-diabatic 
(CD) field [13 . This strategy was recently implemented 
in an atomic Bose-Einstein condensate [ljj. 

In this Letter, we report another experimental imple- 
mentation of the assisted adiabatic passage (AAP), us- 
ing a single nitrogen vacancy (NV)-center in diamond. 
The NV-center system is one of the promising candi- 
date quantum bits (qubits) for building quantum com- 
puters [III, [l]|. Various techniques for implementing 
high-fidelity coherent control of the NV-centers have been 
developed recently [l4l - fl6l |. In our experiment, we im- 
plement the AAP on an electron spin transition of the 
NV-center, using resonant microwave fields as controls. 
The fidelity of the passage is sufficiently high that we 
can perform multiple rounds of the passage in opposite 
directions. The results show good agreement between 



FIG. 1: (color online). Characteristics of the single qubit 
model for A = 0.2 and b = 2. (a) shows the energy levels and 
(b) the time dependence of the three field components. 



theory and experiment. In contrast to the previous im- 
plementation where the coherent control was applied to 
an ensemble of atoms fill], our experiment is implemented 
on a single spin in one NV-center. 

Model.- In close analogy to the Landau-Zener model 

EMU, 

we describe the AAP with the Hamiltonian 



H LZ {t) = \{t)I z + M x , 



(1) 



where the I x x are spin operators, and X(t) and A are 
dimensionlcss fields applied along the z and x direc- 
tions. This Hamiltonian is a model for an arbitrary 
two level system and plays a prominent role in various 
fields of physics, such as in coherent control [U H3] and 
quantum criticality [2lT[23l |. The energy eigenvalues are 
±y/X 2 (t) + A 2 /2, and the instantaneous ground state is 



| fl( i ) ) =sin M|0)_ CO8 M|i) > 



(2) 



where tan#(i) = A/A, 6 E [0, n], and |0) and |1) denote 
the eigenstates of I z with eigenvalues ±1/2, respectively. 
The minimal gap between the two levels is A. Fig. [1] (a) 
illustrates the energy levels for A = 0.2. For the initial 
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FIG. 2: Spectra of the NV electron spin obtained by a) tran- 
sition selective and b) hard pulses. 



discussion, we assume that the control field is scanned 
linearly with time, 



A(t) = &(*-!). 



(3) 



where t € [0,2], and b represents the scan range of A, 
which is [—6, b], and also its rate of change. 

The assisted adiabatic passage can be implemented for 
any non-zero value of A by introducing a CD field that 
is perpendicular to both the x-c omp onent A and the z- 
component A(t) of the field: 0, E3] 



H CD {t) = V CD {t)I v 



where 



v CD {t) 



bA 



[a 2 + \{ty 



(4) 



(5) 



The total field for the AAP is thus a vector B(t) = [A, 
Vco(t), A(f)]. Fig. [U (b) shows the time dependence of 
the three components for b = 2. 

Experimental protocols and results.- For the experi- 
mental test, we used the electron spin of a single NV- 
center in 12 C enriched diamond. The sample used is a di- 
amond single crystal grown at 5.5 GPa and 1400 °C from 
Co-Ti-Cu alloy by using a temperature gradient method. 
As a solid carbon source we used polycrystalline diamond 
plates synthesized by chemical vapor deposition (CVD) 
utilizing 12 C enriched methane. According to the sec- 
ondary ion mass spectroscopy analysis, 12 C enrichment 
of 99.995% was achieved in the grown crystals. The crys- 
tal was irradiated by 2 MeV electrons with a total fluence 
of 10 11 e/cm 2 at room temperature and subsequently an- 
nealed at 1000 °C for 2 hours in vacuum. The reduced 
number of 13 C nuclear spins results in long relaxation 
times, with T 2 * > 100 (is. 

As shown in Fig. [2J the hypcrfine coupling between the 
electron and the 14 N nuclear spin is rj 2.1 MHz. For the 
present experiments, we use the subspace of this system 
that is spanned by the states 7715 = 0, +1 of the electron 
spin, mi = of the nuclear spin. As shown in Fig. [5] a) , 
this subsystem can be excited with excellent selectivity if 
the amplitude of the microwave field remains well below 
the hyperfine coupling constant. We therefore will not 
consider the nuclear spin state in the following. 



For the AAP, the system should be initialized into the 
ground state \g(0)) of -Hlz(O). In the experiment, we ini- 
tialize the system by a laser pulse into the state |0), whose 
overlap with \g(0)) is sin[#(0)/2]. In the experiment, we 
use —3 > A(0)/A > —20, which results in overlaps of 
[0.9871, 0.9997]. For our purpose, this is sufficiently close 
to unity. 

After initialization, the system evolves under the time- 
dependent Hamiltonian 



into the state 



H{t) = H LZ (t) + H CD {t) 



m))=U(t)\g(0)), 



(6) 



(7) 



where U(t) represents the propagator generated by H{t). 
During and after the scan, we read out the state of the 
system by a second laser pulse, which again projects the 
system onto the state |0). We write the probability of 
finding the system in this state as P| ) = |(0|V>(0)| 2 - 

In the actual experiment, the fields A, Vcd and A act 
on the spin in a rotating reference frame. Writing 



= e -m)h 



for the transformation from the laboratory-frame to the 
rotating frame, the Hamiltonians of the two frames are 
related as 



H 1 
H 



lab 



= U r H Uh Ul + iUrU}. 

= UlH TOt U r - iUltj r . 



(8) 



The laboratory-frame Hamiltonian thus has the field 
components 

u x {t) = Acos£(t)+Vc£>(*)sin£(t) 
u y (t) = -Asin£(i) + Vcd (*) cos£ (t) 

dm 



A(f) 



dt 



This Hamiltonian must match the experimentally avail- 
able Hamiltonian, whose general form is 

H cxp = -wo/, + 2wi(t)I x . 

Accordingly, we must have 

e(t) = u t+ f \{t')dt', 

Jo 

which defines our rotating frame transformation. For the 
transverse field components, we invoke the rotating field 
approximation, which allows us to set 

wi(t) = Acos£(i) + V CD {t)sm£{t) 

and ignore the y-componcnt. The amplitude of the field 
is therefore 

Mi) I = VAH^W 2 . 
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Experimental limitations define a maximum possible field 
amplitude, which we designate as f2. In the present ex- 
periment, it is determined by the requirement that no 
transitions of other nuclear spin states are excited, and 
we found a value of 17 ps 2tt ■ 0.2 MHz to be a suitable 
compromise. The maximum field amplitude is reached 
at 

n = ± (9) 

Sa 

where we have defined the scale factor s a , which converts 
the dimensionlcss quantities A and b into actual field 
amplitudes (in Hz) and defines the scan duration 

T S = 2 Sa . (10) 

Figure [3] shows the experimental results of the AAP 
for the scan rates b = 0.6, 1, 1.6, 2, 3 and 4. For these 
parameters, we can approximate s a ~ b/(QA). For b = 
0.6, the scan duration becomes 4.77 /zs and for b = 4 it is 
31.83 [is. The z-component is always scanned from —40 
to +40 kHz (in frequency units), while the y-component 
(the CD field) reaches a maximal amplitude of 200 kHz 
at the anticrossing point. The x-component of the field in 
the rotating frame is 13.3 kHz for b = 0.6 and 2.0 kHz for 
6 = 4. Filled circles show the experimentally measured 
overlaps of the state with the ground state |0). The error 
bars (1 standard deviation) were determined by repeating 
each experiment 10 times. The solid lines, which agree 
very well with the experimental data points, represent 
the theoretical behavior. For comparison, we also show 
one data set that was obtained without CD field, for a 
scan rate of b = 2. These data points are represented by 
the empty circles and the corresponding theoretical curve 
can be approximated as a horizontal line going through 
the experimental points. Clearly, a passage without CD 
field results in an almost completely diabatic transfer. 

A significant additional speedup over this direct imple- 
mentation of the AAP scheme is possible if the amplitude 
of the control field is kept constant throughout the scan, 
at the maximum value that is experimentally accessible. 
For this purpose, we divided the total evoluti on p eriod 
into N segments with identical duration <5 [U, |25| |. The 
propagator U(t) generated by H(t) can then be imple- 
mented by a series of discrete unitary operators U m for 
each period: 

u(js) « n u m . (a) 

m— 1 

For each segment, the propagator U m = e - lSH ( mS ) was 
implemented as a rectangular pulse with the amplitude 
f2, which is the same for all segments. With this ap- 
proach, we expect to increase the possible scan rate. The 
total number N of these segments must be large enough 
that this discrete version becomes a good approximation 
of the continuous scan. 




FIG. 3: (color online). Experimental results obtained with 
the analog implementation of the AAP. The individual ex- 
periments correspond to b = 0.6, 1, 1.6, 2, 3 and 4. The 
experimental data are shown as filled circles, and the corre- 
sponding error bars were obtained by repeating the experi- 
ment 10 times. The curves show the theoretical result for an 
ideal scan. The empty circles with the almost horizontal line 
show the result for a reference experiment without CD field. 

We calculate the Hamiltonians H m for the individual 
segments from the requirement that 

JJ m = e ~ iSH ( mS ) = e -i>T m H m 

where r m is the duration of the pulse that implements the 
m th segment and H m the experimentally realized Hamil- 
tonian in the rotating frame. We can therefore reduce 
the duration of this segment by a factor s m = r m /S 
if we simultaneously scale up the Hamiltonian, H m = 
H(mS)/s m . Using the transformation ([5]), we can calcu- 
late the required laboratory-frame Hamiltonian 

H l £ h = -oj Q Iz + 2CII X cos(uj m t + cj) m ) 

and the required duration r m . The scaling factor 

Sm ~ 5 ~ n [ ' 

is now different for every segment. The angular frequency 
oj m and the phase 4> m become 

AM . , VcpjmS) 
u m = w o 1 , tan0 m = . 

S m A 

Figure |4] shows the experimental results. Here, we used 
the same nominal scan rates b as in the analog case, but 
split the scan into N = 56 segments. The experimental 
data points are represented by filled circles, the error 
bars were obtained by repeating the experiments 8 times. 
Clearly, the experimental data agree very well with the 
theoretical expectation shown as the red curves. The 
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FIG. 4: (color online). Experimental results obtained with 
the rapid-scan approach. For details see the caption of Fig. 
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FIG. 5: (color online). Multiple assisted adiabatic passages, 
(a) shows the time-dependence of the three field components, 
while b) and c) show the experimental data, together with 
the theoretical curves for the ideal case (zero loss). 

empty circles again represent the reference experiment 
obtained by setting Vcd = 0. 

Given the high fidelity obtained in a single passage 
through resonance, we can cycle the system forth and 
back multiple times. The reverse passage is obtained 
by changing b to —6, which changes the direction of the 
scan as well as the sign of the CD field. Figure a) 
shows the time-dependence of all three field components 
in the rotating frame. Figures b) and c) represent the 
experimental results for these repetitive scans, as well 
as the theoretical curves corresponding to the ideal case. 
We find very little loss of population after five passages 
through resonance. 

Discussion. - The experimental implementation of the 
time-dependent Hamiltonian H(t) always occurs with fi- 



nite precision, which results in a loss of fidelity. Exper- 
imental contributions to this loss include the precision 
with which the shaped pulses are implemented - both in 
terms of the amplitude as well as in terms of the time 
resolution. In the rapid-scan approach, the number of 
segments used is an important parameter. We used nu- 
merical simulations of the experiment to estimate these 
losses. The results indicate that for the analog scan, finite 
time resolution of 0.25 ns reduces the fidelity by a frac- 
tion of the order of 10~ 8 . In the rapid-scan approach, for 
N = 56 segments, the maximal loss during the AAP oc- 
curs near the critical point at A = 0. For the lowest scan 
rate, with b = 0.6, this loss is of the order of < 10~ 3 , 
for the faster scan rate, b = 4.0, it rises to 2.8 ■ 10~ 2 . 
However, these are mostly intermediate losses, which are 
recovered during the second part of the scan: The calcu- 
lated loss of fidelity at the end of the evolution period is 
< 10" 4 . 

For the parameters chosen here, the duration of a sin- 
gle scan varies between 4.77 and 31.83 (is in the analog 
version and from 2.0 to 2.5 fxs in the rapid scan imple- 
mentations. They are thus all short compared to the 
coherence time of our sample (T2 ~ 500/is). 

To estimate the speed-up provided by the CD field, we 
used numerical simulations of an unassisted scan, with 
the same parameters as the experimental scan in Fig. 
but slower scan rates. To reach a fidelity of 0.99, the scan 
duration had to be extended to 2.33 ms. This implies 
that the assisted scan allows a speedup of more than two 
orders of magnitude (~ 150) if a linear frequency scan is 
used and of about three orders of magnitude (~ 960) in 
the rapid passage (constant amplitude) version. 

Conclusion.- We have implemented the assisted adia- 
batic passages through analog and rapid-scan approaches 
in a two level quantum system by controlling a single spin 
in a NV-centcr in diamond. This approach allows a sig- 
nificant increase in the scan rate compared to the unas- 
sisted passage and therefore reduces the requirements on 
the decoherence time of the system to which it is applied. 
Like in the unassisted case, the scan has to be slower if 
the minimum gap is small. If the scan is performed lin- 
early in time, the total duration also increases with the 
scan range. However, with the rapid-scan approach that 
we introduced here, the scan range can be increased ar- 
bitrarily with very little time-penalty. Our experiment 
results illustrate the excellent coherent control that can 
be achieved for the spins of NV-centers. These results 
should be helpful for all applications requiring quantum 
adiabatic passages, such as implementing geometric gates 
for quantum computatio n [2 611 , adiabatic control in inter- 
acting two level systems [27j or adiabatic quantum com- 
puting 0| . 

This work is supported by the Hcinrich Hertz Founda- 
tion, and the DFG through grant Su 192/27-1. 
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